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Nicotinic acetylcholine receptor (AChR) channels at neuromuscular synapses rarely open in the absence of ago- 
nists, but many different mutations increase the unliganded gating equilibrium constant (E 0 ) to generate AChRs 
that are active constitutively. We measured E 0 for two different sets of mutant combinations and by extrapolation 
estimated E 0 for wild-type AChRs. The estimates were 7.6 and 7.8 x 10~ 7 in adult-type mouse AChRs ( — 100 mV at 
23°C) . The values are in excellent agreement with one obtained previously by using a completely different method 
(6.5 x 10" , from monoliganded gating). E 0 decreases with depolarization to the same extent as does the dili- 
ganded gating equilibrium constant, e-fold with ^60 mV. We estimate that at —100 mV the intrinsic energy of the 
unliganded gating isomerization is +8.4 kcal/mol (35 kj/mol), and that in the absence of a membrane potential, 
the intrinsic chemical energy of this global conformational change is +9.4 kcal/ mol (39 kj/ mol) . Na + and K + in the 
extracellular solution have no measureable effect on E 0 , which suggests that unliganded gating occurs with only 
water occupying the transmitter binding sites. The results are discussed with regard to the energy changes in recep- 
tor activation and the competitive antagonism of ions in agonist binding. 



INTRODUCTION 

Nicotinic acetylcholine receptors (AChRs) are ligand- 
gated ion channels that regulate the flux of cations across 
cell membranes by switching between closed- and open- 
channel conformations. This global change in AChR 
structure ("gating") occurs both with and without activat- 
ing ligands ("agonists") present at the two extracellular 
transmitter binding sites. However, the gating equilib- 
rium constant increases in the presence of agonists 
because these ligands bind with a higher affinity to the 
open-channel conformation of the protein. 

To understand the energy that agonists provide for 
the full AChR gating isomerization, it is essential to link 
measurements of the equilibrium constants in both unli- 
ganded and fully liganded conditions in a thermody- 
namic cycle (Monod et al., 1965; Karlin, 1967; Jackson, 
1986; Auerbach, 2010). The unliganded gating equilib- 
rium constant (E 0 ) quantifies the tendency of the pro- 
tein to adopt the high affinity/ open conformation 
when no agonists are present at the binding sites, and 
the logarithm of E 0 gives the intrinsic energy of the gat- 
ing conformational change. The extent to which each 
agonist molecule increases the gating equilibrium con- 
stant above this level is a measure of the energy it con- 
tributes toward the global conformational change of 
the protein. Having accurate estimates of this energy is 
one key to understanding how ligand binding and con- 
formational change are coupled. 



In pioneering experiments, Jackson (1984) used single- 
channel electrophysiology of mouse muscle and Neubig 
and Cohen ( 1980) used ion flux measurements of Torpedo 
electroplax membranes to estimate independently that 
the unliganded gating equilibrium constant of these 
AChRs is small, ^10~ 7 . More recently, we estimated E 0 wt 
for adult mouse neuromuscular AChRs (expressed in 
human embryonic kidney [HEK] cells; —100 mV) by 
using two different methods. By measuring single-channel 
currents in receptors having multiple gain-of-function 
background mutations, we estimated E 0 wt ~1.2 x 10~ 7 
at -100 mV (Purohit and Auerbach, 2009). In these 
experiments, energy coupling between mutations was 
not taken into account, so this estimate is low. In the sec- 
ond method, E 0 wt was estimated from the ratio of the 
gating equilibrium constants of diliganded and mono- 
liganded AChRs and was ^6.5 x 10~ 7 at -100 mV (Jha 
and Auerbach, 2010). 

Here, we provide a new and more accurate estimate for 
E 0 wt in adult mouse neuromuscular AChRs (expressed 
in HEK cells), obtained by using an improved version 
of the gain-of-function background approach. We also 
explore the sensitivity of E 0 to membrane voltage and to 
the concentration of monovalent cations in the extra- 
cellular solution. 
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MATERIALS AND METHODS 

Cell culture and mutagenesis 

HEK 293 cells were maintained in Dulbecco's modified Eagle's 
medium supplemented with 10% fetal bovine serum and 1% pen- 
icillin-streptomycin, pH 7.4. The cultures were incubated at 37°C 
and 5% CGv Mutants of mouse AChR subunit cDNAs were made 
by using a site-directed mutagenesis kit (QuikChange; Agilent 
Technologies) and verified by nucleotide sequencing. HEK cells 
were transiently transfected using calcium phosphate precipita- 
tion by incubating them with 3.5-5.5 ug DNA per 35-mm culture 
dish in the ratio of 2:1:1:1 (a/p/8/e) for ~16 h. Most electro- 
physiological experiments were done ^24 h after transfection. 

Electrophysiology 

Single-channel currents were recorded at 23°C in the cell- 
attached patch configuration. The cells were bathed in K + Ringer's 
solution composed of (in mM): 142 KC1, 5.4 NaCl, 1.8 CaCl 2 , 
1.7 MgCl 2 , and 10 HEPES/KOH, pH 7.4. The patch pipettes were 
generally filled with Dulbecco's PBS containing (in mM): 137 
NaCl, 0.9 CaCl 2 , 2.7 KC1, 1.5 KH 2 P0 4 , 0.5 MgCl 2 , and 8.1 Na 2 HP0 4 , 
pH 7.3/NaOH. In experiments to study the effect of ions (Figs. 6 
and S4), the pipette solution contained (in mM): 0-150 NaCl/ 
KCl/CsCl, 0.9 CaCl 2 , 0.5 MgCl 2 , 1.5 KH 2 P0 4 , and 8.1 Na 2 HP0 4 . 
To study the effect of near complete ion removal on the receptor, 
the patch pipettes were filled with H 2 0 and 0.1 mM CaCl 2 . In 
some experiments, the partial agonist choline (Cho) was added 
to the patch pipette solution. 

Patch pipettes were fabricated from borosilicate glass, coated 
with sylgard (Corning), and fire polished to a resistance of ^10 
Mft when filled with pipette solution. Single-channel currents 
were recorded using an amplifier (PC505; Warner Instruments). 
The single-channel records were low-pass filtered at 20 kHz 
and digitized at a sampling frequency of 50 kHz using a data 
acquisition board (SCB-68; National Instruments). The pipette 
holder used for the unliganded experiments was never exposed 
to agonists. 

Selection and characterization of background 

In wild-type (wt) AChRs, unliganded openings are very rare, so 
Eq" 1 cannot be measured directly. We engineered two different 
sets of gain-of-function background mutation combinations to 
study unliganded gating. The first set of mutants, most of which 
had the background combination DYS (aD97A + aY127F + 
ctS269I), were measured and described in detail previously 
(Purohit and Auerbach, 2009). In the second set, the key muta- 
tion was at txA96. Mutations here can increase the diliganded 
gating equilibrium constant with Cho (E 2 cho ) substantially and 
over a wide range (Cadugan and Auerbach, 2010). In addition, 
we used 0T456I, 8I43Q eE181T, eL269F, and 8V269A as back- 
ground mutations in four different combinations. Individually, 
each of these mutations increases E 2 cho (Table SI) : (3T456I, 2.1-fold 
(Mitra et al., 2004); 8I43Q 5.0-fold (see supplemental text); 
eE181T, 2.2-fold (Jha et al., 2012); eL269F, 179-fold (Jha et al., 
2007); and 8V269A, 250-fold (Cymes et al., 2002), in all cases by 
an almost equivalent increase in E 0 . 

We assumed that the effects of the individual mutations, as 
measured by E 2 cho , were caused only by equivalent fold changes in 
E 0 and were approximately energetically independent of each 
other. For example, the product of the individual fold increases 
in E 2 cho for the background (PT456I + 8I43Q+ eL269F + eE181T) 
is 4135, and the individual fold increases for the oA96L and 
F mutants are 50 and 497. With the above assumptions, adding 
these to the background should result in a total increase in E 0 by 
factors of ~0.21 and 2.1 million. The E 0 obs values for these final 
constructs were 0.11 and 1.37, which are approximately in the 



same ratio. This indicates that the energetic consequences of the 
c«A96 mutations were approximately independent of those of the 
background. 

Expression of most of the aA96 mutants was robust. In the ab- 
sence of any agonists, the currents were ^7 pA at —100 mV, and 
there was no channel block. 

Kinetic modeling and determination of rate constants 

Kinetic analyses of single-channel currents were done by using 
the QuB software suite (http://www.qub.buffalo.edu). For esti- 
mation of rate constants, clusters of individual channel activity, 
flanked by >20-ms nonconducting periods, were selected by eye. 
Clusters of unliganded activity were idealized into noise-free 
intervals after they were low-pass filtered digitally (12 kHz) by 
using the segmental K-means algorithm (Qin et al., 1996). The 
unliganded opening (f 0 ) and closing (bo) rate constants were 
estimated from the idealized intra-cluster interval durations by 
using the maximum-interval likelihood (MIL) algorithm after 
imposing a dead time of 50 us (Qin et al., 1997). 

The kinetics of unliganded gating activity are complex 
(Jackson et al, 1990; Grosman and Auerbach, 2000b; Purohit and 
Auerbach, 2009) . There were multiple components in the intra- 
cluster interval duration histograms. The idealized intra-cluster 
intervals were first fitted by a simple two-state closed (C)-open 
(O) model using the MIL algorithm. Additional C and O states 
were added, one at a time, until the log-likelihood score failed 
to improve by >10 units. The unliganded opening and closing 
rate constants (f 0 and b 0 ) were estimated from the inverse life- 
times of the predominant (>80%) component of closed- and 
open-interval duration histograms. E 0 , the unliganded gating 
equilibrium constant, was calculated as the ratio fo/b 0 . The fact 
that the slope of the E 0 plot (Fig. 2) was equal to 1 indicates that 
the above method consistently identified the same unliganded 
gating rate constants. 

The unliganded activity of some oA96 background combina- 
tions exhibited distinct modes. In such cases, clusters were sepa- 
rated into populations based on their mean open probability, 
mean closed time, and mean open time using the "Select" fea- 
ture and the K-means algorithm of QuB. The predominant 
(^70%) gating mode was selected and E 0 was determined as 
described above. 

Voltage dependence analysis 

For unliganded gating, E 0 was measured at different membrane 
potentials (Fig. 5, A-D). Two different methods were used to 
estimate the voltage dependence of E 2 . In some experiments, 
we assumed that low affinity agonist equilibrium dissociation 
constant was independent of voltage and measured an apparent 
diliganded opening rate constant (f 2 *) and the diliganded clos- 
ing rate constant (b 2 ) using 50-100 uM Cho. At these concentra- 
tions, channel block by Cho was minimal. In other experiments, 
we estimated f 2 using 20 mM Cho and b 2 using 100 uM Cho 
(Fig. 5, G and H). 

To study the perturbation of voltage on gating, we plotted log f 0 
(or f 2 ) versus log E 0 (or E 2 ) for eight different constructs (Fig. 5). 
This is a rate-equilibrium (R-E) relationship. The slope of this 
plot (<t>) was estimated by an unweighted linear fit in Sigmaplot 
10 (Systat Software, Inc.). For a linear, sequential, bounded reac- 
tion mechanism, <1> may provide relative temporal information 
regarding the movement of the perturbed side chain in the gat- 
ing reaction ( 1 is "early," 0 is "late," and the same is "synchronous") 
(Zhou etal., 2005). 

Online supplemental material 

The online supplemental material contains Tables S1-S5 and 
Figs. S1-S4. Tables SI and S2 illustrate the effect of mutant com- 
binations on the gating equilibrium constant. Tables S3 and S4 
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show the voltage dependence of the unliganded (E 0 ) and dili- 
ganded (E 2 *) gating rates and equilibrium constants. Table S5 deals 
with the effect of cations on E 0 . Fig. SI shows the locations of the 
mutated amino acids in the Torpedo AChR (Protein Data Bank 
accession no. 2bg9) . Fig. S2 presents the high resolution view of 
example unliganded currents for different aA96 mutants. The 
effects of membrane voltage and extracellular Cs + ions are shown 
in Figs. S3 and S4. The online supplemental material is available 
at http://www.jgp.org/cgi/content/full/jgp.201110752/DCl. 



RESULTS 

Theory and assumptions 

A cyclic reaction scheme provides a useful, if approxi- 
mate, description of AChR activation (Fig. 1). With- 
out an external energy source, the product of the 
equilibrium constants linking the unliganded closed 
channel (R) and the diliganded open channel (A 2 R*) 
is the same for the clockwise and anticlockwise con- 
necting paths. Hence, with two equivalent binding 
sites (Akk et al., 1996; Salamone et al., 1999; Jha and 
Auerbach, 2010), 



E 0 X 



(1) 



where E 2 is the diliganded gating equilibrium constant; 
E 0 is the unliganded gating equilibrium constant; and X 
is the square root of their ratio or, equivalently, the 
ratio of low and high affinity equilibrium dissociation 
constants, Kd/J d . The logarithm of X. quantifies the en- 
ergy provided by each agonist molecule to power gating. 

Many mutations throughout the AChR increase E 2 , as 
measured using the partial agonist Cho (Grosman and 
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Figure 1 . Cyclic reaction scheme for AChR binding and gating. 
A is the agonist, R is the resting conformation of the AChR (low 
affinity, closed channel), and R* is the active conformation (high 
affinity, open channel). Rj and J d are the equilibrium dissociation 
constants of R and R*. E 0 , E 1( and E 2 are the gating equilibrium 
constants with zero, one, and two bound agonist molecules. The 
logarithm of E 0 is the intrinsic energy for the protein isomeriza- 
tion in the absence of bound agonists. In wt AChRs, the unboxed 
states are rarely visited and only the boxed states determine physi- 
ological and pharmacological responses. The two transmitter 
binding sites are assumed to be functionally equivalent. Without 
an external energy source, the product of the equilibrium con- 
stants for the outer paths connecting R with A 2 R* is equal (E 2 /Kd 2 = 
Eo/Jd 2 ) , so E 2 = E 0 A. 2 where X = Kj/Jd. A mutation that changes E 2 
can do so by changing E 0 , X, or both. 



Auerbach, 2000a; Chakrapani et al., 2004; Mitra et al., 
2004; Spitzmaul et al., 2004; Purohit et al., 2007). 
According to Eq. 1, a mutation can increase E 2 by in- 
creasing E 0 , X, or both. The evidence suggests that many 
mutations that are not in the immediate vicinity of the 
transmitter binding sites that increase E2 Cho do so by an 
approximately equivalent increase in E 0 and have little 
or no effect on X (Jadey et al., 2011). These same experi- 
ments also show that these mutations often influence E 0 
independently, because for many different mutant 
combinations, the net fold change in E 2 cho is approxi- 
mately the product of the fold changes for each muta- 
tion alone. 

These two characteristics of many AChR mutations — 
that they only change the intrinsic gating energy and 
act independently — offer a simple way to estimate the 
intrinsic gating energy of wt AChRs. First, measure the 
fold increase in E 2 cho in AChRs having a side-chain sub- 
stitution that increases the frequency of unliganded 
openings and assume that the increase is caused exclu- 
sively to a parallel increase in E 0 (that the mutation did 
not alter X cho ). Second, calculate what the predicted 
fold increase in E 2 would be for a combination of muta- 
tions (E 2 mut /E 2 wt ) assuming that the effect of each muta- 
tion is energetically independent. That is, assume that 
the fold change in E 0 for the mutant combination is the 
product of the fold changes for each mutation alone. 
Third, express AChRs having the mutant combination 
and measure experimentally the unliganded gating 
equilibrium constant, E 0 obs . If the two assumptions are 
perfectly correct, then E 0 obs /E 0 wt = E 2 mut /E 2 wt . A log-log 
plot of E 0 obs versus E 2 mut /E 2 wt should be linear with a 
slope of 1 and an intercept on the ordinate axis, where 
log (E 2 mut /E 2 wt ) = 0, which is the estimate of log (E 0 wt ). 

Estimates of Eo™* 

We used two different sets of background mutations to 
estimate E 0 wt in adult-type AChRs (Fig. SI). The first 
set was the same as reported previously but analyzed 
without considering the possibility that there was energy 
coupling between the mutations (Purohit and Auerbach, 
2009) . Many of the mutated amino acids in the first set 
were in the a subunit, close to each other and the 
transmitter binding site. For the second set, we chose 
fewer and more separated amino acids. The key muta- 
tion of the second set was at residue aA96. Side-chain 
substitutions here increase E 2 cho over a wide range 
and, in some cases, to a great extent (Cadugan and 
Auerbach, 2010). As a consequence, only a few mutations 
were needed to generate the large increase in unli- 
ganded gating required to produce the clusters of sponta- 
neous openings that are necessary for measuring E 0 obs . 

Fig. 2 shows single-channel records of unliganded 
gating activity from example mutants of the second set. 
Without any agonists in the pipette or bath, openings 
occurred in clusters that each reflect the constitutive 
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gating activity of a single AChR. The gaps between clus- 
ters reflect epochs when all of the AChRs in the patch 
were desensitized. Different aA96 mutation increased 
the unliganded open probability of clusters to different 
extents. E 0 obs was measured for each mutant combina- 
tion from the durations of shut and open intervals 
within clusters (see Materials and methods) . 

Log-log plots of E 0 obs versus E 2 mut /E 2 wt for both mu- 
tant sets are shown in Fig. 2 B. In both cases, the results 
were well described by a straight line over approximately 
four orders of magnitude. For the first mutant set, the 
extrapolated value was E 0 wt = 7.8 x 1(T 7 (95% confidence 
interval, 1.6-37 x 10~ ), and for the second set, it was 
E 0 wt = 7.6 x 1(T 7 (3.1-18 x 1(T 7 ). The slope for the first 
set was 0.85 ±0.1, and for the second it was 0.99 + 0.07. 

The slope for the first set was <1, which indicates 
that the two core assumptions of the method were not 
perfectly satisfied. In the first set, there was a small 
amount of energy coupling between the mutations, a 
small effect of some mutations on \ cho or both. Because 
the previous analysis of these mutants did not take 
into account such coupling, the new E 0 wt estimate was 
^6.5 times higher than the previous value. The slope 
of correlation for the second set was ~1, so we con- 
clude that there was essentially no energy coupling 
between these mutations. 



Fig. 3 A shows unliganded gating activity from adult- 
type AChRs having only the aA96H substitution in 
both a subunits. This side-chain substitution alone was 
sufficient to raise unliganded gating activity to a level 
where clusters of spontaneous openings from individ- 
ual AChRs were apparent. Fig. 3 B shows that adding 
a distant mutation (in the M4 helix of the (3 subunit) 
increased E 0 obs further and, hence, the probability of 
being open within unliganded clusters. This mutation 
increased E 2 cho by increasing the opening rate constant 
by ^T. 1-fold and decreasing the closing rate constant by 
^2. 1-fold (Mitra et al., 2004). In our experiments with 
unliganded AChRs, we found a similar result, namely 
that this mutation increased E 0 obs by increasing the 
opening rate constant by 1.8-fold and decreasing the 
closing rate constant by 2.2-fold. 

In the absence of agonists, the intra-cluster interval 
duration distributions show multiple closed and open 
components (Jackson, 1986; Grosman and Auerbach, 
2000b; Purohit and Auerbach, 2009). Recently, these 
"extra" states have been incorporated into a gating 
model that posits that the transmitter binding sites 
are energized ("primed") before channel opening, with 
brief versus long unliganded openings reflecting AChRs 
having one versus both sites so activated (Mukhtasimova 
etal., 2009). It is therefore important to consider whether 



A txA96V (no agonist, +(3T456I+sE181T+eL269F) 



2 s 



n 



8 1 in i r TUT 



IT 



B 





3 




2 




1 




0 




-1 


<n 
.a 
o 

O 


-2 


UJ 






-3 


o 






-4 




-5 




-6 




-7 



aA96L 



500 ms 



aA96V 




01 23456789 10 11 
expected log (E 2 mut /E 2 wt ) 

352 Unliganded gating of AChRs 



Figure 2. Estimating E 0 M using mutations 
that increase the gating equilibrium constant. 
There was no agonist present in the pipette 
or bath. (A) Low time-resolution view of clus- 
ters of unliganded openings (open is down; 
— 100 mV). Each cluster reflects the constitu- 
tive gating activity of an individual AChR. The 
background was (JT456I + eE181T + eL269F 
(Table SI). (B) Estimation of the unliganded 
gating equilibrium constant of wt AChRs 
(Eo™ 1 ) from different mutant sets (1, open; 2, 
closed). The data from set 1 are from Puro- 
hit and Auerbach (2009). The x-axis value is 
the log of the product of the fold increases in 
E 2 cho for each mutation alone (Table SI), and 
the y axis is the log of E 0 ° bs measured directly 
from the clusters. For both mutant sets, the 
extrapolated ordinal intercept was approxi- 
mately —6.12 (unweighted straight-line fit 
shown with 95% confidence intervals), which 
corresponds to E 0 wt = ^7.6 x 10~ 7 . (Inset) 
Higher time-resolution view of clusters from 
mutant set 2 (aA96) expressed on the same 
background. 



or not our association of the main component of each 
unliganded interval duration distribution with sojourns 
in the unliganded states of the cyclic model (R<-»R* in 
Fig. 1) was appropriate. 

Two previous results bear on this assessment. First, 
many point mutations of aromatic residues at the trans- 
mitter binding site eliminate the extra components ap- 
parent with the first mutant set and reduce unliganded 
gating to a simple two-state mechanism (Purohit and 
Auerbach, 2010). For example, the mutation ctY198F, 
which removes only two oxygen atoms from the ^300-kD 
AChR complex and has essentially no effect on diliganded 
gating, results in unambiguous R<-»R* unliganded gat- 
ing. Indeed, 55 of 57 mutations of the a-subunit-binding 
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Figure 3. "Extra" shut and open components in unliganded gat- 
ing. In all panels, there was no agonist present in the pipette or 
bath. (A; top) Continuous trace of clustered unliganded openings 
from AChRs having only the aA96H mutation (in both a subunits) . 
This single mutation increases substantially E 0 and, hence, constitu- 
tive activity. (Bottom) Interval duration histograms of events within 
clusters and an example cluster (underlined, above). Two com- 
ponents are apparent in both the shut and open histograms. E 0 obs 
was estimated from the predominant components. (B) Adding the 
[5-subunit M4 mutation (3T456I increases Eq° s but does not elimi- 
nate the extra components. (C) Adding the binding site mutation 
c«W149M eliminates the long openings but has little effect on E 0 obs . 
(D) Another example of extra components in both the shut- and 
open-interval duration histograms. (E) These extra components 
disappear with the addition of the binding site mutation oiW149F. 
After correcting for the backgrounds, the E 0 values were the same 
with or without the extra components. 



site aromatic residues (oY93, aY190, aY198, and aW149) , 
and so far only of these four amino acids, completely 
eliminate the long component of unliganded openings. 

We observed similar results using the ctA96 back- 
ground constructs. Fig. 3 C shows that with a construct 
from the second mutant set, the extra open component 
disappears when the oW149M mutation is added to the 
background. However, the E 0 obs estimate was similar with- 
out (0.04) or with (0.02) this added mutation. Fig. 3 D 
shows that multiple shut and open components are 
clearly apparent in the oA96Y background (used in 
mutant set 2) having the added (3M4 mutation, but that 
these disappear with the addition of the aW149F muta- 
tion. Again, when the effect of the aW149F on E 0 is taken 
into account, the E 0 obs values were the same without and 
with (0.02) this mutation. Hence, for both the first and 
second mutant sets, the extra shut and open unliganded 
components are reduced by mutations of oW149, with- 
out influencing E 0 obs estimates. 

A second previous observation regarding the extra 
unliganded components is that they are absent in cur- 
rents generated by liganded AChRs (see Fig. 4 A) . E 0 wt was 
estimated by examining ligand-gated currents from 
AChRs that had only one functional transmitter binding 
site (Jha and Auerbach, 2010). The results from these 
experiments, where no extra components were apparent 
and gating activity was unambiguously AR<->AR*, yielded 
a similar estimate of E 0 wt (for two different agonists) to 
that obtained from unliganded AChRs (6.5 x 10~ 7 ). 

Several additional results support the hypothesis that 
the main components of unliganded gating reflect the 
same conformational change as in diliganded gating 
and, hence, can be linked in a thermodynamic cycle. 
The liganded and unliganded R* conductances are the 
same (Jackson, 1984). The voltage dependences are the 
same (Fig. 5) . Gating values of many different residues 
are approximately the same with versus without ligands, 
which suggests that the conformational pathway connect- 
ing the R and R* end states is approximately the same 
(Purohit and Auerbach, 2009) . The energy estimates cal- 
culated by using the cycle for the two ACh molecules 
( — 10.2 kcal/mol) as well as the intrinsic isomerization 
(+8.4 kcal/mol; see below) are physically realistic. Collec- 
tively, the observations indicate that the association of the 
main shut and open components of unliganded gating 
with the R<->R* states of the cyclic model is appropriate. 

Voltage dependence of the gating equilibrium constant 

We next measured the voltage dependence of the gat- 
ing equilibrium constant in both diliganded and unli- 
ganded AChRs (Fig. 4) . Although the voltage dependence 
of the open-channel lifetime and the channel-closing rate 
constant in diliganded AChRs have been measured 
many times previously, estimates of the voltage depen- 
dence of the diliganded equilibrium constant are rarer 
(Colquhoun and Sakmann, 1985; Grosman etal., 2000) . 
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To measure both the forward and backward gating 
rate constants in unliganded AChRs, we often needed 
to use background mutations that increased the open- 
ing rate constant and decreased the closing rate con- 
stant. In combination, these made it possible to measure 
intra-cluster closed- and open-interval durations over 
the range of - 120 to +60 mV. 

The results of these experiments are shown in Fig. 5 
(Table S3). For both unliganded and diliganded condi- 
tions, and regardless of the background mutations, the 
gating equilibrium constant changed exponentially 
with membrane potential, decreasing on average e-fold 
with an ^59 ± 6-mV depolarization (mean + SD) . These 
observations indicate that the background perturba- 
tions we used did not alter significandy the voltage sen- 
sitivity of the receptor. We estimate that in the absence 
of a membrane potential, E 0 wt ' 0 mV = 1.3 x 10~ 7 . 

We were also able to measure the extent to which the 
membrane potential changes the gating equilibrium 
constant by altering the forward versus the backward 
rate constant. The insets in Fig. 5 show this analysis in 
the form of R-E plots. For both the diliganded and unli- 
ganded conditions, and for all backgrounds, the pre- 
dominant effect of depolarization was to increase the 
channel-closing rate constant. The average value of the 
slope of the eight R-E plots (<t>) for the voltage sensor 
was 0.34 ±0.11. 

Ion dependence of E 0 

In the above experiments, the "unliganded" condition 
was AChRs exposed to the pipette solution without ago- 
nists, which was PBS (that contained 137 mM NaCl). 



Many AChR agonists are cations, so we considered the 
possibility that some or all of the constitutive activity ob- 
served in PBS might be attributed to extracellular Na + 
occupying the binding sites and acting as an agonist 
(X. > 1 ) . To test this hypothesis, we removed all mono- 
valent cations from the pipette solution (which con- 
tained only water and 0.1 mM CaCy and clamped the 
cell membrane potential to +70 mV so that the currents 
flowed through the channel in an outward direction. 
The AChRs had the mutations oA96Y (to increase the 
frequency of unliganded openings) and 8V269A (to 
prolong the open-channel lifetime) . Fig. 6 A shows that 
AChRs still gave rise to clusters of spontaneous openings 
in the monovalent cation-free solution. The E 0 of the 
mutant in this condition was not significantly different 
from the E 0 observed under normal conditions with 
137 mM Na + (Table S5). We exposed the AChRs in the 
patch to different concentrations of monovalent cations 
(Na + , K + , or Cs + ) in the pipette solution (the bath was K + 
Ringer's solution; see Materials and methods). Fig. 6 B 
shows that the pipette [Na + ] had no effect on E 0 , and that 
the same result was observed when the Na + was replaced 
with K + . The effect of pipette Cs + is described in the 
online supplemental material (Fig. S4 and Table S5). 

DISCUSSION 

All of the estimates of E 0 wt , obtained with different 
methods and different sets of background mutations, 
are in excellent agreement. Given the error limits on 
these measurements, we think a reasonable (and easily 
remembered) value for adult-type mouse neuromuscular 
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Figure 4. Interval duration histograms and example currents at different membrane potentials. (A) Diliganded gating activated by 
50 pM Cho (the construct was ctA96L + 8V269A) . (B) Unliganded gating (no agonist in the pipette or bath; the construct was aA96Y + 
8V269A) . The extra components in the interval duration histograms in unliganded gating are not apparent in diliganded gating. Both 
with and without agonists, the shut lifetime decreases and the open lifetime increases with hyperpolarization. 
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AChRs is E 0 wt = 7 x 10~ 7 ( - 100 mV at 23°C) . Mutations 
throughout the AChR can change E 0 , so this value 
cannot be applied directly to any other receptor, for 
example, from a different species or one comprised of 
different subunits. 

The logarithm of an equilibrium constant is propor- 
tional to the free energy difference between the end 
states. From the logarithm of E 0 wt , we estimate that at 
— 100 mV, the intrinsic energy for the unliganded gating 
isomerization is +8.4 kcal/mol (35 kj/mol). The electric 
field of the membrane contributes a modicum of energy 
for this conformational change. The gating equilibrium 
constant decreases e-fold with an ^60-mV depolariza- 
tion, so at zero membrane potential, E 0 wt = 1.3 x 10~ 7 . 
We estimate that the intrinsic chemical energy of the 
mouse adult-type neuromuscular AChR-unliganded gat- 
ing isomerization is +9.4 kcal/mol (39 kj/mol). 

With the new E 0 wt estimate, we can update the key equi- 
librium constants for adult AChR activation at —100 mV 



(Fig. 1). E 2 ACh « 25 (Wang et al., 2000; Hatton et al., 
2003; Jadey et al., 2011), so, from the relationship A. = V 
(E 2 /E 0 ) , we estimate that X ACh ~ 6,000. Each ACh mol- 
ecule increases the gating equilibrium constant by this 
factor and provides approximately —5.1 kcal/mol 
(21.4 kj) to power the gating isomerization of the pro- 
tein. Together, the two binding sites provide —10.2 
kcal/mol, which is sufficient to overcome the +8.4- 
kcal/mol intrinsic energy difference and, hence, open 
the channel with a high probability. With regard to 
transmitter binding, the low affinity equilibrium dis- 
sociation constant is K/ 011 ~ 166 uM (Jadey et al., 201 1) , 
so we calculate (from the relationship A. = K<i/J d ) that 
the high affinity constant is J d ACh = 28 nM. Most of the 
extra energy from the agonist in the open conforma- 
tion of the binding site is in the form of heat (Gupta 
and Auerbach, 201 1 ) . 

We can separate the intrinsic gating free energy of 
unliganded gating AG E0 = +8.4 kcal/mol into enthalpy 
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Figure 5. Voltage dependence of the gating 
equilibrium constant. (A-D) The unliganded 
gating equilibrium constant E 0 as a function of 
the membrane voltage (V m ) for different con- 
structs. (Inset) 5> value analysis of the voltage 
sensor. (E-H) The diliganded gating equilib- 
rium constant E 2 as a function of V m for different 
mutants. Overall, there was an e-fold decrease in 
the gating equilibrium constant with a depolar- 
ization of 60 mV with <1> = 0.34. Unliganded con- 
structs (no agonist): A, aDY + eL269F + eP245L; 

B, aDYS + 8L265T and aDY + eL269F + eP245L; 

C, aA96Y + (3T456I + 8I43Q + eE181T + eL269F; 

D, oA96Y + 8V269A. Diliganded constructs (acti- 
vated by Cho): E, oA96L + 8V269A; F, eE181T + 
eL269F; G, eL269F; H, eS450A. 
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(AH) and entropy (AS) components. The temperature 
dependence of E 2 ACh estimates the sum of the enthalpy 
changes for both the transmitter molecules and the 
intrinsic isomerization: AH E2(ACh) = AH ACh + AH E0 . The 
experimental value was AH E2(ACh) = +0.7 kcal/mol 
(Gupta and Auerbach, 201 1) . These same experiments 
showed that for different agonists, the differences in AG 
for diliganded gating was essentially equal to the differ- 
ences in AH. Therefore, we postulate that AH ACh = - 10.2 
kcal/mol. We now calculate that AH E0 = +10.9 kcal/mol, 
and, from the relationship AG = AH — TAS and the value 
AG E0 = +8.4 kcal/mol, that TAS E0 = -2.5 kcal/mol. 

Many AChR mutations alter E 0 without changing 
or \, so many mutation-induced changes in cellular 
dose-response and kinetic profiles, including some that 
cause congenital myasthenic syndromes, can be traced 
specifically to changes in E 0 (Zhou et al., 1999). Indeed, 
it is remarkable that a single mutation in loop A of the 
otsubunit (aA96H) increases E 0 by ^1 00,000-fold (-6.8 
kcal/ mol) and by itself can lead to a relatively high fre- 
quency of opening in the absence of activating ligands 
(Fig. 3). This result is consistent with the previous sug- 
gestion that this region of loop A is acting as a "latch" 
that serves to hold the AChR in the low affinity, closed- 
channel conformation (Chakrapani etal., 2004). 

The voltage dependence of E 0 is approximately the 
same as for E 2 (Fig. 5). The gating equilibrium constant 
decreases e-fold with an ^60-mV depolarization, with or 
without agonists present at the binding sites. This result 
indicates that the product of the charge and displace- 
ment of the AChR voltage sensor is the same regardless 
of whether or not ligands occupy the transmitter binding 
sites. For all of the mutant combinations, the €> value of 
the voltage sensor was low and similar with or without 



agonists, which is consistent with the proposal that the 
sequence of conformational changes within the com- 
plete gating isomerization is similar in liganded versus 
unliganded AChRs (Purohit and Auerbach, 2009). 

We found no evidence that extracellular Na + or K + 
influences E 0 . When added to the extracellular solu- 
tion, these cations were without effect on unliganded 
gating (Fig. 6) . Indeed, E 0 was the same without any 
monovalent cations as it was in PBS. We hypothesize 
that constitutive gating of AChRs occurs with only water 
present at the transmitter binding sites. 

Extracellular Na + and K + are, however, competitive in- 
hibitors of ACh association (Akk and Auerbach, 1996) . 
There are two ways to reconcile this result with the lack 
of effect of these ions on E 0 . It is possible that Na + and 
K + do indeed occupy the transmitter binding sites but 
experience no change in binding energy when the pro- 
tein changes between its low and high affinity forms; 
i.e., perhaps \ Na+ = 1. This seems improbable. Other 
charged ligands experience a large energy change when 
the binding site switches its conformation. For example, 
tetramethylammonium ion experiences an approxi- 
mately — 4.0-kcal/mol energy change (X TMA = 879) in 
this switch (Zhang et al., 1995; Jadey et al., 2011). Also, 
the resolution of our experiments is approximately ±0.5 
kcal/mol, so \ Na+ must be almost exacdy 1 (between 0.4 
and 2.3) to have no measureable effect on E 0 , which 
is unlikely. Regardless, the above conclusion that unli- 
ganded gating occurs in water only is correct even if Na + 
(and K + ) does occupy the binding sites but has X ~ 1. 

A second possibility is that the competitive inhibition 
between ion and agonist occurs outside the binding site, 
before the low affinity complex has formed. Recently, we 
have found that low affinity binding to the AChR requires 
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Figure 6. Effect of Na + and K + on E 0 . (A) Sponta- 
neous currents (no agonist in pipette or bath) at 
+70 mV (open is up). The pipette solution was only 
water plus 0.1 mM Ca 2+ . (Bottom right) Marked 
cluster shown at higher resolution. (B) Effect of 
adding monovalent cations to the pipette solution. 
(Left) Unliganded current clusters with different 
amounts of extracellular Na + . (C) The cluster open 
probability does not change with [Na + ]. (Right) E 0 
was not influenced by [Na + ] or [K + ] . 
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both diffusion and a protein conformational change 
(unpublished data) . This implies that there is an "encoun- 
ter complex" between agonist and protein that exists after 
the ligand has diffused to the binding site but before the 
conformational change has taken place. It is possible that 
Na + competes with ACh in the formation of this encoun- 
ter complex rather than in low affinity binding. Further 
experiments are needed to reveal the basis of the compe- 
tition between extracellular cations and agonists. 
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